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(Dated: January 8, 2009) 

A search for pair production of the lightest supersymmetric partner of the top quark, ti, is per- 
formed in the lepton+jets channel using 0.9 fb _1 of data collected by the DO experiment. Kinematic 
differences between titi and the dominant top quark pair production background are used to sepa- 
rate the two processes. First limits from Run II of the Fermilab Tevatron Collider for the scalar top 
quark decaying to a chargino and a 6 quark (ti — > Xib) are obtained for scalar top quark masses of 
130-190 GeV and chargino masses of 90-150 GeV. 
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Supersymmetry [l| introduces a superpartner for each 
of the left and the right-handed top quarks. Because of 
the large top quark mass, the mixing between those two 
can be substantial and lead to a large difference in the 
mass eigenvalues of the two scalar top ("stop") quarks. 
Thus, the lighter stop quark t\ could possibly be the 
lightest scalar quark and within reach at the Fermilab 
Tevatron Collider. In the Minimal Supcrsymmetric Stan- 
dard Model (MSSM) stop quarks are produced mainly in 
pairs \t\t\) via the strong interaction, the same mech- 
anism as for top quark pair production (ti) [2J. The 
expected next-to-leading-order (NLO) cross section at 
a center of mass energy of 1.96 TeV for a stop quark 
of mass equal to 175 GeV is (0.58±g;^) pb [3], while 
for a top quark of the same mass the cross section is 
(6.8±0.6) pb [4|. The stop quark pair production cross 
section strongly depends on the mass of the stop quark. 

The different possible decay modes of the stop quark 
result in a number of distinct final state signatures. The 
branching ratios for stop quark decays depend on the 
parameters of the model, in particular the masses of the 
supersymmetric particles involved. The decays to a c 
quark and the lightest neutralino (ti c\i) [5j and to a 
b quark, a lepton, and a sneutrino {t\ — > bfc^iif) Q have 
already been explored at DO in Run II of the Tevatron. 
For stop quarks lighter than the top quark the decay 
channel t\ — > x^fo, with subsequent decay of the lightest 
chargino xt to tne lightest neutralino Xi an d a W boson, 
can dominate, if kinematically allowed. In this Letter we 
assume that the branching ratio B(t\ — > Xi&) = 1. This 
channel has been explored only once before by the CDF 
collaboration in Run I of the Tevatron at y/s = 1.8 TeV 
for stop quark masses of 100-120 GeV ■ With a dataset 
more than ten times larger, we obtain first limits in this 
channel at s/s = 1.96 TeV for stop quark masses in the 
range 130-190 GeV. 

The titi event signature in the studied decay chan- 
nel can be similar to the ti signature, making it possi- 
ble for the titi signal to be embedded in the ti event 
sample. The goal of this analysis is to search for this 
possible hidden admixture. The main difference relative 
to ti production is the additional presence of neutralinos 
in the event. However, this does not lead to significantly 
higher missing transverse energy ($r)> since the neutrali- 
nos are mostly produced back-to-back. We consider the 
decay channel with one W boson decaying to hadrons 
and the other one to an electron or muon and a neutrino. 
Scenarios with both on-shell and off-shell W bosons pro- 
vide the same signature. The final state consists of one 
high-pr lepton, Iflx from the neutrino and the neutrali- 
nos, two jets originating from b quarks ( "6 jets" ), and two 
light-quark jets. This is referred to as the "lepton+jets" 
channel. We consider twelve mass points, for which the 
studied decay can dominate. We fix the neutralino mass 
to 50 GeV, a value close to the experimental limit from 
LEP ||, and we vary the stop quark mass from 130 to 



190 GeV and the chargino mass from 90 to 150 GeV to 
obtain the desired event signature. For larger neutralino 
masses the signature changes and the sensitivity of this 
study decreases. 

The search is conducted using data collected by the 
DO detector [9J in pp collisions at the Fermilab Tevatron 
Collider. Triggers require an electron or muon and at 
least one jet with large transverse momentum (px)- The 
dataset comprises an integrated luminosity of 913 pb -1 
for events containing electrons in the final state, and 
871 pb -1 for events with muons. 

We select events with one isolated electron with px > 
20 GeV and pseudorapidity \rj\ < 1.1, or one isolated 
muon with p T > 20 GeV and \r]\ < 2.0, and fx > 
20(25) GeV in the electron (muon) channel [1(3]. To re- 
ject events with mismeasured leptons, the lepton momen- 
tum vector and the Jpx vector are required to be sepa- 
rated in azimuth. In addition, we only accept events 
with >3 jets, each with px > 15 GeV and \rj\ < 2.5, 
of which the jet with largest px ("leading jet") has to 
have pr > 40 GeV. Events with a second isolated elec- 
tron or muon with px > 15 GeV are rejected. Details 
about object identification, jet energy corrections, and 
trigger requirements can be found in Ref. [nj. In addi- 
tion, we require at least one 6-tagged jet in each event, 
where the b jets are identified through a neural network 
algorithm 11 1. 

For events with four or more jets, a kinematic fitting 
algorithm is used to reconstruct the objects to a it 
hypothesis, which is used to separate i±ti from tt events. 
The fitter minimizes a x 2 statistic within the constraints 
that both candidate W boson masses are 80.4 GeV and 
that the masses of the two objects reconstructed as top 
quarks are the same. The fitter considers only the four 
jets of highest px, uses 6-tagging information to minimize 
combinatorics, and varies the four-vectors of the detected 
objects within their resolution. Only events for which 
the fit converges (86-95% of signal events depending on 
the mass point and lepton flavor) are selected for further 
analysis. 

The events are classified into four distinct subsamples, 
according to jet multiplicity (3 jets or > 4 jets) and lepton 
flavor (e+jets or /i+jets). All subsamples are used to 
obtain the final limit. 

Because of their topological similarity to the signal, ti 
events are the most challenging background. Of the other 
background processes, production of W bosons in associ- 
ation with jets (W+jets), and multijet events, where jets 
are misidentified as isolated leptons, are most important. 
Far smaller contributions arise from Z+jets, single top 
quark, and diboson production. 

Except for the multijet background, the shape of dis- 
tributions in all processes are modeled through Monte 
Carlo (MC) simulation. The t\t\ signal is generated by 
pythia v6.323 [IH in its general MSSM mode, where the 
top trilinear coupling A t and the SU(2) gaugino mass M2 
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TABLE I: Expected numbers of events with total uncertain- 
ties and observed numbers of events after the final selection. 



Sample 


=3 jets 




> 4 jets 






e+jets 


/i+jets 


e+jets 


/x+jets 


Signal 

m h [GeV]/m 


-±[GeV] 








190/150 
130/90 


q n+0.3 
°- z -0.3 

io.4±i;2 


9 9+0.2 
z - z -0.2 

6 5+ 06 

°-°-0.8 




9 1+0.3 
n 9+O.6 



Background 

ti 

W+jets 
Z+jets 
Single top 
Diboson 
Multijet 
Total 



1-16.0 
-15.3 
«7 7+25.5 
Ol./ _24.4 
1.5 
1.1 
1.6 
1.5 
1.1 
0.9 

22.3±4.2 

-,+31.4 
-29.1 



77.6Z 



5.2^ 
9.3^ 
4.2 4 



186.2Z 



ro r + 12.0 
JO-J_19.6 

TTA±%\ 
6.91 2 ; 3 
7 5+ 13 

1.2 

q Q+ 1.0 

J - 8 - 0.9 

3.0±2.4 
1 57 2+ 22 ' 9 

10 r -^_29.6 



103.0tgj 



1-12.1 
-12.8 
5+ 0.8 
0.7 
|+ 0.7 
0.7 
0.5 
0.3 

10.7±2.6 

1+26.8 
-26.6 



17.1] 
2.1 

3.V 
1.4" 1 



138.1] 



84.211?;° 
21.61 7 7 . 8 o 

q q+ 0.9 

°- J - 0.8 
9 C+ 0.7 
z -°- 0.6 
1 9+ 0.4 
0.3 

3.3±2.7 

11fin +20.9 
IIP. U_ 20. 9 



Data 



193 



163 
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are varied to set the stop quark mass and the chargino 
mass, respectively. The neutralino mass is kept at the 
same value by keeping the U(l) gaugino mass Mi con- 
stant. The ti background is also generated by PYTHIA, 
using a top quark mass of 175 GeV. The VF+jets and 
■Z+jets processes are generated by ALPGEN 2.05 [lj] for 
the matrix element calculation, with subsequent par- 
ton showering and hadronization generated with PYTHIA. 
Single top quark events are generated by CompHEP- 
Singletop [IH and diboson production is modeled by 
PYTHIA. All generated events are passed through a 
GEANT-based 16] simulation of the DO detector and re- 
constructed using the same software as for data. To im- 
prove agreement between data and MC simulation, addi- 
tional corrections [Io| are applied to the simulation before 
selection. 

The contribution of the multijet background for each 
jet multiplicity and lepton flavor is determined from data 
using a method which exploits the fact that this back- 
ground contains jets that mimic leptons, whereas the 
other processes have a true isolated lepton [T3| ■ The nor- 
malization of the T^+jets background is estimated be- 
fore imposing the 6-tagging requirement, by subtracting 
from data: (i) the estimated multijet background, and 
(ii) the ti, Z+jets, single top, and diboson contributions 
as calculated from their next-to-leading order cross sec- 
tions [3, fla |. The remaining events are assumed to be 
Vy+jets background, where we have scaled the heavy fla- 
vor component (Wbb plus Wcc) by a relative factor of 
1.17±0.18. This factor was derived on a statistically in- 
dependent sample with two jets and at least one 6-tag. 

Table U shows the numbers of expected and observed 
events after the final selection, found to be in good agree- 
ment. For signal events the mass points with the highest 
and lowest event yield are shown as examples. 

Because of the similarity of the iiti and ti final 
states [191, a multivariate likelihood discriminant |20l is 



employed to discriminate between the two processes. We 
study the kinematic differences and choose the variables 
of greatest discrimination and low correlation as input 
to the multivariate discriminant. For events with three 
jets, where the two jets besides the leading 6-tagged jet 
are referred to as light jets j, the following five vari- 
ables are used: (i) the invariant mass of the three jets, 
(ii) K^ n = AR$ n pf in , where Ai?™ in is the minimum 
AR [lfj separation between a pair of jets (in rapidity- 
azimuth space) and p™ m is the minimum jet px in that 
pair, (iii) the smaller of the AR separations between the 
leading 6-tagged jet and either the lepton or the vector 
sum of the two light jets, (iv) the px of the system of the 
two light jets, and (v) the lepton- J$t transverse mass [21 1. 
For events with four or more jets, the following five vari- 
ables are used: (i) the top quark mass as reconstructed 
by the kinematic fitter, (ii) the scalar sum of the pt of 
the four leading jets, (iii) the invariant mass of the sys- 
tem of the second and third leading jet, excluding the 
leading 6-tagged jet, (iv) if™ ln , and (v) the pr of the 
fourth leading jet. 

Figure [T] shows the variable with the greatest separa- 
tion for each jet multiplicity as a comparison between 
data and the prediction. Figure [2] shows the resulting 
discriminant for the mass point with 771^=175 GeV and 
ni~± — 135 GeV in the 3-jet and the 4-jet subsample, com- 

paring the prediction with data. The prediction for t\t\ 
signal (solid line) peaks at 1, while it peaks at for ti. 

We use a Bayesian approach [22j | to extract upper 

limits on the stop quark pair production cross section 

(a r r ) from the discriminant distributions. We construct 
*i*i 

a binned likelihood as a product over all bins in the dis- 
criminant distribution as well as each of the four chan- 
nels considered, assuming a Poisson distribution for the 
observed counts per bin. For the signal cross section, 
we assume a flat non-negative prior probability. By in- 
tegrating over signal acceptance, background yields and 
integrated luminosity using Gaussian priors for each sys- 
tematic uncertainty, we obtain the posterior probability 
density as a function of cross section for signal. The up- 
per limit on <7j j at 95% confidence level is the point 
where the integral over the posterior probability density 
reaches 95% of its total. 

We differentiate between systematic uncertainties that 
change the yield uniformly for all bins of the discrim- 
inant, and those that affect each bin differently. The 
effects are given as a percentage on the event yield of the 
affected process; they can vary widely, depending on the 
subsample and the physics process. The sources chang- 
ing the yield uniformly include the uncertainties on in- 
tegrated luminosity (6.1%) [23}, efficiency of the event- 
based data quality selection (0.5%), theoretical cross sec- 
tions (13-20%), top quark mass (1.3-7%), estimation of 
the VF+jets background (24-74%, depending on the jet 
multiplicity and lepton flavor subsample), influence of 
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FIG. 1: Comparison of the prediction with data after the final selection for the e+jets and /x+jets channels combined, a) the 
invariant mass of the three jets in events with 3 jets, b) the reconstructed top quark mass in events with >4 jets. The solid 
line shows the distribution for a signal point, enhanced by a factor of ten. 
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FIG. 2: Comparison of the discriminant distribution for data with the prediction after the final selection for the e+jets and 
^t+jets channels combined, for events with a) 3 jets, and b) >4 jets. The solid line shows the distribution for a signal point, 
enhanced by a factor of ten. 



the signal on the W+jets normalization (0.8-3.4%), es- 
timation of the multijet background (19-84%, depend- 
ing on the subsample), lepton identification and recon- 
struction efficiencies (2.2-2.5%), primary vertex identi- 
fication efficiency (2.7%), and trigger efficiencies (1.2- 
2.7%). The sources that also change the shape of the 
discriminant distribution include jet energy scale calibra- 
tion (0.6-30%), and 6-tagging (0.1-27%). Limits on the 
stop quark pair production cross section are degraded by 
about a factor of two when all systematic uncertainties 
are accounted for. 

Table [II] shows the results for each mass point for the 
combination of all channels. The results are also illus- 
trated in Fig. [3] The expected limits are derived from 
the sum of all selected background samples without a t±ti 
contribution, but including the tt background according 
to its theoretical cross section. The observed limits on 



the cross section are a factor of 2 — 13 larger than the the- 
ory prediction and agree with the expected limits within 
uncertainties. In some cases, most notably for the mass 
point with mi — 175 GeV and m~± — 135 GeV, the ob- 

1 Al 

served limit is higher than the expected limit, pointing 
to an excess of signal-like data. To quantify the signif- 
icance, the peak position of the posterior probability is 
compared to its width. In this case, the peak is 1.62 
standard deviations away from zero. 

In summary, we present first limits on the t\ti pro- 
duction at the Tevatron Run II for a light stop quark 
of 130-190 GeV decaying to a b quark and the lightest 
chargino. In the MSSM scenarios studied by this search, 
we derive upper limits on the cross section that are a fac- 
tor of 2 — 13 above the theory prediction and agree with 
the expected limits within uncertainties. 

We thank the staffs at Fermilab and collaborating 
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TABLE II: The predicted titi cross section and the expected 
and observed Bayesian upper limits on the titi cross section 
at the 95% confidence level for different assumed values of 
and m~±. We assume m ; o = 50 GeV and B(ti — > yfb) = 1. 
The uncertainties on the theoretical prediction result from the 
simultaneous variation by a factor of two of the factorization 
and renormalization scales about their nominal values, set 
equal to the stop quark mass. The uncertainties on the ex- 
pected limits represent the one standard deviations estimated 
via background-only pseudo-experiments. 
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